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1. Introduction 

The use of a generalized oscillator algebra for characterizing a dynamical system gave rise to a 
great deal of papers. Among many works, we may quote the polynomial Heisenberg algebra worked 
out in the context of supersymmetry [1-3], the deformed Heisenberg algebra introduced in connection 
with parafermionic and parabosonic systems [4-7], the CVextended oscillator algebra developed in the 
framework of parasupersymmetric quantum mechanics [8-11], and the generalized Weyl-Heisenberg 
algebra Wk related to Z^. -graded supersymmetric quantum mechanics [12-16]. In this direction, the 
construction of a truncated generalized oscillator algebra was developed by several authors. In particular, 
the pioneer work along this line by Pegg and Barnett led to calculating the phase properties of the 
electromagnetic field [17]. Let us also mention the works [18,19] in relation with orthogonal polynomials 
of a discrete variable and [16] in connection with phase operators and dynamical systems. 

Recently, a generalized oscillator algebra A K , a one-parameter algebra that is a particular case of 
the algebra W\, was studied for the purpose of defining phase operators and the corresponding phase 
eigenstates [16]. In addition, it was shown that the phase states for A K with k > 0, which are particular 
coherent states [20,21], can serve to construct mutually unbiased bases which are of considerable interest 
in quantum information and quantum computing [16]. 

It is the aim of the present paper to analyze the algebra A K from the point of view of group 
theory. Since A K can describe the Morse system for k < as well as the harmonic oscillator and 
the Poschl-Teller systems for k > 0, we expect that the groups SU(2) and SU(1,1) play a central role. 
The search for phase operators and temporally stable phase states thus amounts to study generalized 
coherent states for SU(2) and SU(1,1). 

The material presented here is organized as follows. Section 2 deals with the generalized oscillator 
algebra A K and its connection with the Lie algebra of SU(2) and SU(1,1). The phase operators and 
the phase states introduced in [16] are described in the framework of SU(2) and SU(1,1). Section 4 is 
devoted to a truncation of the algebra A K . In section 5, the phase operator for the group SU(2) is shown 
to be of relevance for the determination of mutually unbiased bases (cf. [22-32]). Finally, the quadratic 
transformation that connects the phase states for SU(2) to angular momentum states is studied in 
Section 6. This transformation generalizes the discrete Fourier transform whose the main properties 
are given in the appendix. 

The notations are standard. Let us simply mention that: 5 a> b stands for the Kronecker symbol of a 
and b, I for the identity operator, for the adjoint of the operator A, and [A, B] for the commutator 
of the operators A and B. The bar indicates complex conjugation and matrices are generally written 
with bold-face letters (I d is the <i-dimensional identity matrix). We use a notation of type \ip) for a 
vector in an Hilbert space and we denote and \4>){^\ respectively the inner and outer products of 

the vectors \tp) and \<f>). As usual N, N*, Z and M + are the sets of integers, strictly positive integers, 
relative integers and positive real numbers; K. and C the real and complex fields; and Z/rfZ the ring of 
integers 0, 1, . . . , d — 1 modulo d. 
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2. Generalized Oscillator Algebra 

2.1. The Algebra 

Following [16], we start from the algebra A K spanned by the three linear operators a", a + and N 
satisfying the following commutation relations 



a , a 



I + 2kN, [N,a ± }=±a ± , (a~) f = a+ = N (1) 



where k is a real parameter. In the particular case k = 0, the algebra *4. is the usual harmonic oscillator 
algebra. In the case k / 0, the operators a - , a + and iV in (1) generalize the annihilation, creation 
and number operators used for the harmonic oscillator. Thus, the algebra A K can be referred to as 
a generalized oscillator algebra. In fact, the algebra A K represents a particular case of the generalized 
Weyl-Heisenberg algebra Wu introduced in [12-15] to describe a fractional supersymmetric oscillator. A 
similar algebra, namely the Ca -extended oscillator algebra, was studied in connection with a generalized 
oscillator [8-11]. 

2.2. The Oscillator Algebra as a Lie Algebra 

The case k = corresponds of course to the usual Weyl-Heisenberg algebra. It can be shown that the 
cases k < and k > considered in [16] are associated with the Lie algebras of the groups SU(2) and 
SU(1,1), respectively. We shall consider in turn the cases when k < and k > 0. 

For k < 0, we introduce the operators J_, J + and J3 defined via 

J- := -^a~, J + := — , J 3 := — (I + 2kN) (2) 



2k 



-K 

They satisfy the commutation relations 

[J + ,J_]=2J 3 , [J 3 ,J+] = J + , [J 3 ,J_] = -J_ (3) 

and therefore span the Lie algebra of SU(2). 

Similarly for k > the operators K-, K + and K 3 , given by 

K_:=-^a~, K + :=-^a + , K 3 := -L (/ + 2kN) (4) 

\/K \/K 2K 



lead to the Lie brackets 

[K+, K-] = -2K 3 , [K 3 , K+] = K + , [K 3 , K_] = -K. (5) 
of the group SU( 1,1). 

2.3. Rotated Shift Operators for Su(2) and Su( 1,1 ) 

We are now in a position to reconsider some of the results of [16] in terms of the Lie algebras 
su(2) and su(l,l). This will shed new light on the usual treatments of the representation theory of 
SU(2) and SU(1,1) as far as the action on the representation space of the shift operators of these groups 
are concerned. 
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Let us first recall that in the generic case (k e M), the algebra A K admits a Hilbertian representation 
for which the operators a~, a + and N act on a Hilbert space T K spanned by the basis {\n) : n = 0, 1, . . .} 
that is orthonormal with respect to an inner product (n\n') = 8 n>n >. The dimension of T K is finite when 
k < or infinite when k > 0. The representation is defined through [16] 

a + |n) = ^/F(n + l)e~ i[F{n+1) - F{n)]ip \n + 1} (6) 
ar\n) = ^/F^n)e +i[F{n) - F{n - l)]v \n - 1) (7) 
a"|0)=0, N\n)=n\n) (8) 

where if is an arbitrary real parameter and the function F : N — > M + satisfies 

F{n + 1) -F(n) = 1 + 2ku, F(0) = ^ F(n) = n[l + k{u - 1)] (9) 

Obviously, for k > the dimension of JF K is infinite. In contrast, for k < the space is 
finite-dimensional with a dimension given by 

d := 1 - - with - - e W (10) 

K K 

It is thus possible to transcribe (6)-(8) in terms of the Lie algebras su(2) and su(l,l). 
2.3.1. The Su(2) Case 

Let us consider the (2j + l)-dimensional irreducible representation of SU(2) spanned by the 
orthonormal set 

B 2 j+i ■= {\j,m) :m = j,j-l,..., -j} (11) 
where \j, m) is an eigenvector of J 3 and of the Casimir operator 

J 2 := J+J- + J 3 (J 3 -1) (12) 

We know that 

J 2 \j,m) =j(j + l)\j,m), J 3 \j,m) = m\j,m) (13) 
with m = j, j — 1, . . . , — j for fixed j (2j e N). Following [25-32], we make the identifications 

|n) •f-)- \j, m), n -h- j + m (14) 



Consequently, we have 



which leads to the relation 



d = 2j + \ = \-- (15) 

K 



2jk = -1 & ~- = 2j (16) 

K 
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that is crucial for the connection between A K and su(2). It is to be noted that (14)— (16) are compatible 
with (2). We can then rewrite (6) and (7) in the su(2) framework. In fact, by combining (2), (9) and (16) 
with (6) and (7), we obtain 



J+\j,m) = y/(j -m)(j + m + l)e~ 2im ^\j, m + 1) (17) 
J-\j,m) = + m)(j - m + l)e 2i(m - 1)KV |j, m - 1} (18) 

Equations (17) and (18) differ from the usual relations, well known in angular momentum theory, by the 
introduction of the phase factor (p. The standard relations, that correspond to the Condon-Shortley phase 
convention of atomic spectroscopy, are recovered when ip = 0. 

Although there is no interdiction to have ip ^ 0, it is worthwhile to look for the significance of 
the introduction of <p. Let us call J + and J_ those operators J + and J_ which correspond to ip = 0, 
respectively. It is easy to show then that J± and J± are connected by the similarity transformation 

J ± = e ~ lXKip J±e lXKV (19) 

where the operator X reads 

X := J 2 - J 3 (J 3 - 1) = J + J_ = J + J_ (20) 

Note that the nonlinear transformation J± J±, defined by (19), leaves invariant the Casimir operator 
J 2 of SU(2). We shall see in section 5 that the parameter <p is essential in order to generate mutually 
unbiased bases. 

2.3.2. The Su( 1,1) Case 

The representation theory of SU(1,1) is well known (see for example [20]). We shall be concerned 
here with the positive discrete series D' + of SU(1,1). The representation associated with the Bargmann 
index k can be defined via 



K+\k, k + n) = ^(2k + n)(n + l)e~ i ^ k ' n) \k, k + n + l) (21) 

K_\k, k + n) = ^{2k + n-l)ne^ n - 1] \k, k + n — 1) (22) 

K 3 \k,k + n) = (k + n)\k,k + n) (23) 

with 

K 2 \k,k + n) = k(l - k)\k,k + n), K 2 := K + K_ - K 3 (K 3 - 1) (24) 

where n E NandfC 2 stands for the Casimir operator of SU( 1,1). This infinite-dimensional representation 
is spanned by the orfhonormal set {\k, k + n) : n E N}. Equations (21) and (22) differ from the standard 
relations [20] by the introduction of the real-valued phase function ip. Such a function is introduced, 
in a way paralleling the introduction of the phase factors in (6) and (7), to make precise the connection 
between A K and su(l,l) for n > 0. The relative phases in (21) and (22) are such that K + is the adjoint 
of K-. For fixed k and k, we make the identification 

|n) <+\k,k + n) (25) 
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Then, from (4) we get the central relation 

2kn = 1 <=> - = 2k (26) 

K 

to be compared with (16). Furthermore, by combining (4), (6), (7), (21), (22), (25) and (26) we get 

Fin) = n 



1+ 2^- i; 



${h,n) = T {h + n)(p (27) 
k 



Finally, the action of the shift operators K + and K_ on a generic vector \k,k + n) can be rewritten as 

K + \k, k + n) = y/(2k + n){n + l)e- 2i{k+n)K ^\k, k + n + l) (28) 
K_\k, k + n) = v /(2k + n-l)ne 2i{k+n - 1)KV \k, k + n - 1} (29) 

The particular case (p = in (28) and (29) gives back the standard relations for SU(1,1). 

The operators K + and K_ are connected to the operators K + and K- corresponding to <p = by 

K± = e lYnv K±e-' iYKip (30) 

with 

Y := K 2 + K 3 (K 3 — 1) = K + K_ = K + K_ (31) 

so that the nonlinear transformation K± -H- K±, defined by (30), leaves invariant the Casimir operator 
K 2 ofSU(l,l). 

3. Phase Operators 

Phase operators were defined in [16] from a factorization of the annihilation operator or of A K . We 
shall transcribe this factorization in terms of the lowering generators J_ and K_ of SU(2) and SU(1,1), 
respectively. 

3.1. The Su(2) Case 

Let us define Ed via 

J_ = E d yfj^JZ (32) 

The operator E d can be developed as 

E d = ** imr - 1)mp \j t m-l)(j,m\ (33) 

m=-j 

where m — 1 should be understood as j when m = —j. Consequently 

E d \j, m) = e 21 ^" 1 )^^', m - 1) for m^-j (34) 

and 

E d \j, -j) = e~^\jj) for m = -j (35) 
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It is clear from (34) and (35) that the operator E d is unitary. 

In order to show that E d is a phase operator, we consider the eigenvalue equation 



j 



E d \z) = z\z), \z) := J2 d m zi +m \j,m), z E C, d m G C (36) 



m=-j 



It can be shown that the determination of normalized eigenstates \z) satisfying (36) requires that the 
condition 

z 2j+1 = 1 (37) 
be fulfilled. Hence, the complex variable z is a root of unity given by 

z = q a , q = e 2 ^ +1 \ a = 0,l,...,2j (38) 

As a result, the states \z) depend on a continuous parameter (p and a discrete parameter a. They shall be 
written as \(p,a). A lengthy calculation leads to 

1 j 

\z) = \<p,a) = . e i{j+m)U - m+1)KV q a{j+m) \j,m) (39) 

V J m= _j 

The latter states satisfy 

E d \<p,a) =q a \<p,a) = e 2n[a/ ^ +1) \^a), a = 0,l,...,2j (40) 

Thus, the states \ip,at) are phase states and the unitary operator E d is a phase operator, with a 
non-degenerate spectrum, associated with SU(2). Furthermore, the eigenvectors of E d satisfy 

U(t)\<p,a) = \<p + t,a) (41) 

where 

U(t) := e~' lH \ H := -kX = -kJ + J_ (42) 

and t is a real parameter. Equation (41) indicates that the phase states \(p,a) are temporally stable, an 
important property to determine the so-called mutually unbiased bases [16]. Note that they are not all 
orthogonal (the states with the same y? are of course orthogonal) and they satisfy the closure property 

^\<p,a)(<p,a\ = I (43) 



a=0 



for fixed if (see also [16]). 
3.2. The Su(l,l) Case 
By writing 



K_ = E^^K^Kl (44) 
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it can be shown that 

oo 

= Y,e 2i{k+n)Kip \k,k + n)(k,k + n + l\ (45) 



n=0 

The operator E^ has the following property 



^oo^oo) 1 " = (^oo) f ^oo + | A;, k)(k, k\ = I (46) 



Thus, it is not unitary in contrast with the case of the operator Ed for su(2). 
Let us look for normalized states \z) such that 



oc 



E 00 \z)=z\z), \z) :=^c n z n \k } k + n), zGC, c„6C (47) 



n=0 

One readily finds that 



r) = y/l - \z\ 2 Y,z n e~' m{2k+n ~ 1)Kip \k, k + n), \z\ < 1 (48) 



n=0 

up to a phase factor. Following [33] and [16], we define the states \ip, 9) by 



where 9 G [— it, +7r[. One thus obtains that 



V9,^):= lim t====== 1z) (49) 



1^,5)=^ e i^ e -in(2fc+n-l)^| A , ; fc + ^ (50) 
n=0 

The states (50), defined on the unit circle S 1 , have the property 

E^ip^} =e w \ip 1 9} 1 -tt<9<ti (51) 

The operator E^ is thus a nonunitary phase operator associated with SU(1,1). As a particular case of the 
phase states \(p, 9), the states |0, 9) corresponding to tp = are identical to the phase states introduced in 
[33] for SU(1,1). The parameter (p ensures that the states \(p, 9) are temporally stable with respect to 

U(t) := e~ m \ H := kY = kK + K_ (52) 

in the sense that 

u{t)\ip,e) = \<p + t,e) (53) 

for any real value of t. Note that, for fixed cp, the phase states \ip, 9), satisfy the closure relation 

— / d9\ip,9)((p,6\ = I (54) 



2tt . 

but they are neither normalized nor orthogonal 
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4. Truncated Generalized Oscillator Algebra 

The idea of a truncated algebra for the harmonic oscillator goes back to Pegg and Barnett [17]. 
Truncated algebras for generalized oscillators were introduced in [16,18,19]. In [16], a truncated 
oscillator algebra A KjS associated with the algebra A K was considered both in the infinite-dimensional 
case (k > 0) and the finite-dimensional case (k < 0). The introduction of such a truncated algebra makes 
it possible to define a unitary phase operator for k > and to avoid degeneracy problems for k < 0. We 
shall briefly revisit in this section the truncation of the generalized oscillator algebra A K in an approach 
that renders more precise the relationship between A K>S and A K . 

Let us start with the two operators 

c+ = a + -^^Fi^e-' 1 ^-^- 1 ^^)^-^ (55) 

n=s 

d(K) 

c~ = a--^v^)e +i[i?(n) - F(n - 1)l¥, |n-l)(n| (56) 

n=s 

where d(K) = d — 1 or oo according to whether n < or k > 0. The finite truncation index s is arbitrary 
for k > and less than d for k < 0. It is straightforward to prove that 

c + \n) = ^fF{^+Y)e-' l[F{n+l) - F{n)] *\n + 1) for n = 0, 1, . . . , s - 2 (57) 

c + \n) = for n — s — 1, s, . . . , d{n) (58) 

c~\n) = ^/F{^)e +i[F[n) - F ^- 1)]Lp \n-l) for n = 1, 2, . . . , s - 1 (59) 

c~\n) =0 for n = and n = s, s + 1, . . . , d(n) (60) 

Therefore, the operators c~ and c + = (c~)t lead to the null vector when acting on the vectors of the 
space T K that do not belong to its subspace T R>S spanned by the set {|0), |1), . . . , \s — 1)}. In this sense, 
c + and c~ differ from the operators b + and b~ of [16]. 

In the light of Equations (57)-(60), the passage from the algebra A K to the truncated algebra A KiS 
should be understood as the restriction of the space T K to its subspace JF K a together with the replacement 
of the commutation relations in (1) by 

d(K) 

[ c -,c + ]= I + 2kN - F(s)\s-l)(s-l\-J2(l + 2Kn)\n){n\, [N,c ± }=±c ± (61) 

n=s 

which easily follow from (55) and (56). It should be observed that the difference between the operators 
and b^ 1 manifests itself in (61) by the summation from n = s to n = d(n). 

5. Mutually Unbiased Bases 

5.1. Quantization of the Phase Parameter 

We now examine the consequence of a discretization of the parameter y? in the su(2) case (k < 0). By 
taking (cf. [16]) 



Symmetry 2010, 2 1470 
the state vector \ip, a) becomes 

1 j 

lip, a) = \aa) = , V gtf+^-™+W2+(;+m)«u m \ (63) 

y/2j + 1 ^— ' . 
v m=—j 

The phase operator is of course (^-dependent. For the quantized values of (p given by (62), 
Equations (34) and (35) can be rewritten as 

E d \j, m) = q {m - 1)a \j, m - 1) for m ^ -j (64) 

and 

E d \j, -j) = q ja \j,j) for m = -j (65) 

The corresponding operator E d is thus a-dependent. However, the eigenvalues of E d do not depend on a 
as shown by (40). 

5.2. Connecting the Phase Operator with a Quantization Scheme 

The eigenvector \aa) of E d given by (63) is a particular case, corresponding to r = 0, of the vector 

i j 

\ja;ra) = = ^ q U+rn)U-rr+l)a/2-j™+V+m)a^ m } (66) 

v m=—j 

obtained from a polar decomposition of su(2) [22-26]. More precisely 

\aa) = \ja; 0a) (67) 

In quantum information, \aa) can represent a qudit in dimension d = 2j + 1. The case of a qubit 
corresponds to d = 2, to an angular momentum j = 1/2. 
The vector \ ja\ ra) is an eigenvector of the operator 

i-i 

"ra := e 2 ^ r \j,-j)(j,j\ + q {j - m)a \j,m+l)(j,m\ (68) 

m=-j 

where reR and a G Z/(2j + 1)Z. The action of v ra on | j, m) reads 

tValj, m> = d„,,< 2:: > Jr ./• -j) + (1 - 5 mj )g°" m)a b', m + 1) (69) 
and the matrix elements of v ra in the basis B 2 j+i are 

(j, m\v ra \j, rri) = >),„, /„/y 2;:! ' r ./. -j) + £ m ,m'+i(l - 5 mld )q( j - m >\j, m + 1) (70) 

where m,m' — j,j — 1, . . . , — j. 

As a matter of fact, we have the eigenvalue equation 

v ra \ja; ra) = q^ +a) - a \ja; ra), a = 0, 1, ... , 2j (71) 
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The spectrum of v ra is not degenerate. The vectors \ ja\ ra) are common eigenvectors of J 2 and v ra . For 
fixed r and a, they satisfy the orthogonality relation 

(ja;ra\j/3;ra) = S a> p (72) 

for a, f3 — 0, 1, ... , 2j. 

The operator f ra is unitary and it commutes with the Casimir operator J 2 of SU(2). The set {J 2 , v ra ] 
is a complete set of commuting operators that provides an alternative to the scheme {J 2 , J z }, used in 
angular momentum theory. In other words, for fixed j, r and a, the set 

B ra :={\ja;ra): a = 0,1,..., 2j} (73) 

constitutes a nonstandard orthonormal basis for the (2j + 1) -dimensional irreducible representation of 
SU(2). The basis B ra is an alternative to the canonical basis B 2 j + i defined in (11). The reader may 
consult [22,23] for a study of the {J 2 , v ra ] scheme and of its associated Wigner-Racah algebra. 

The a-dependent operator Ed and the operator v ra are closely connected. Indeed, it can be 
checked that 

E d = q> a (v 0a y = e 2 ^ 2 ^ 1 ) (v 0a y (74) 

as can be guessed from (40) and (71). 

5.3. Introduction of Mutually Unbiased Bases 

The case r = deserves a special attention. Let us examine the inner product (aa\b(3) of the vectors 
| a/3) and \b/3) defined by (63), in view of its importance in the study of mutually unbiased bases (MUBs). 
For a = b, we have 

{aa\af3) = 5 at/ 3 (75) 

Therefore, for fixed j and a (2j e N and a in the ring TLj (2j + 1)Z), the basis 

B 0a :={\aa) : a = 0,1, . . . ,2j} (76) 

(a particular case of the basis B ra ) and the basis B 2 j+i are interrelated via 

(j,m\aa) = — 1 q U+m)(j-m+i)a/2+( j+m )a \(j,m\aa)\ = , 1 (77) 

u 1 1 yjlj + 1 lv 1 71 V2JTT 

with a — 0, 1, . . . , 2j and m — j, j — 1, . . . , — j. In view of (77), we see that B 0a (and more generally 
B ra ) can be considered as a generalized Fourier transform of B 2 j+i. 

For a b, the inner product (aa\b/3) can be expressed in term of the generalized quadratic Gauss sum 
defined by (see [34]) 

\w\-l 

S{u,v,w) := e i7r{uk2+vk)/w (78) 

fc=0 
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In fact, we have 

(aa\bfi) = j^S(u,v,w) (79) 

where 

u:=a-b, v :=-(a-b)d-2(a- P), w := d = 2j + 1 (80) 

The sum S^w, u, w) can be calculated in the situation where u, v and w are integers such that u and w 
are mutually prime, uw is not zero, and uw + v is even. 

Let us now briefly discuss the reason why (63) is of interest for the determination of MUBs. We recall 
that two orthonormal bases of the rf-dimensional Hilbert space C d are said to be unbiased if the modulus 
of the inner product of any vector of one basis with any vector of the other one is equal to 1/ y/d [35,36]. 
For fixed d, it is known that the number N MUB of MUBs is such that 3 < N MUB < d + 1 and that the 
limit N MUB = d + 1 is attained when d is a power of a prime number [35,36]. Then, equation (77) shows 
that any basis B 0a (a G Z/ (2j + 1)Z) is unbiased with B 2 j + i for arbitrary value of 2j + 1. Furthermore, 
in the special case where 2j + 1 is a prime integer, the calculation of S(u, v, w) with (80) leads to 

| M ^)| = __L_ (81) 

for a b, a = 0, 1, . . . , 2j and (3 = 0,1,..., 2j. Equation (81) implies that B 0a and B b for a and b in 
the Galois field F 2:/ +i are mutually unbiased. 

Thus one arrives at the following conclusion. For 2j + 1 prime, the 2j + 1 bases B 0a (a = 0, 1, ... , 2j) 
and the basis B 2 j+i form a complete set of d + 1 = 2j + 2 MUBs. This result is in agreement with the 
one derived in [24-32]. It can be extended to the case r ^ as follows. For arbitrarily fixed r and 2j + 1 
prime, the 2j + 1 bases B ra (a = 0, 1, ... , 2j) and the basis B 2 j+i form a complete set of d + 1 = 2j + 2 
MUBs. The parameter r serves to differentiate various families (or complete sets) of MUBs. 

6. Discrete Fourier Transforms 

We discuss in this section two quadratic versions of the discrete Fourier transform (DFT), namely, 
the quantum DFT that connects state vectors in an Hilbert space and the classical DFT used in signal 
analysis. 

6.1. Quantum Quadratic Discrete Fourier Transform 

Equation (66) shows that the vector \ ja; ra) can be considered as a quantum DFT that is quadratic (in 
m) for a ^ 0. This transform is nothing but a quantum ordinary DFT for r = a = [37]. For fixed j, r 
and a, the inverse transform is 

1 23 

\j, m) = q -U+rn)U-rn+i)a/2+j m r___Y^ q-^ +m ^ a \ja; ra) (82) 

Compact relations, more adapted to the Fourier transform formalism, can be obtained by going back to 
the change of notation given by (14) and (15). Then, Equations (66) and (82) read 

\ja; ra) = qi^^jL y q n(d-n)a/2+n[a-(d-l)r/2] ^ (g3) 



Symmetry 2010, 2 



1473 



and 



d-l 

\ n ) = q -n(d-n)a/2-(d-lfr/4+n(d-l)r/2j_ g-«*|j a; ra ) (g4) 



We shall put 



(F ra ) nQ := J_g"( rf -™)«/2+(<i-l) 2 r/4+n[a-(d-l)r/2] ^ 

\fd 



or 



(F ra ) nQ = -^ e 2 " i//d with / := ±(d - l) 2 r + ±[2a + da - (d - l)r\n - \an 2 (86) 
Vd 

a relation that defines (for fixed d, r and a) ad x d matrix F ra . Let us recall that for a fixed value of d in 
N*, both r and a have a fixed value (reK and a G Z/dZ) and n, a = 0, 1, . . . , d — 1. 
For d = 2j + 1 arbitrary, we can show that 

((F ra ) t F sb ) a ^ = (ja; ra|j/3; sb) (87) 

Therefore, in the particular case r = s and d = p, where p is prime, we have 

\((F ra yF rb ) a p\ = {(j^raljfcrtyl = \(aa\b/3)\ = for a ^ b (88) 

Equation (88) shall be discussed below in terms of Hadamard matrices. 

6.2. Quadratic Discrete Fourier Transform 

6.2.1. Factorization of the Quadratic DFT 

We are now prepared for discussing the transforms (83) and (84) in the language of classical signal 
theory. Let us consider the transformation 

i = {i w GC:m = 0,l d-l} <-> y = {y n G C : n = 0, 1, . . . , d - 1} (89) 

defined by 

d-l d-l 

f« = ( F ra) m „^m ^ X m = ^ ( F ra) mn 2/n (90) 

m=0 n=0 

The particular case r = a = corresponds to the ordinary DFT. For a ^ 0, the bijective transformation 
x -H- y can be thought of as a quadratic DFT. The analog of the Parseval-Plancherel theorem for the 
ordinary DFT can be expressed in the following way. The quadratic transformations x <H- y and x' -H- y' 
associated with the same matrix F ra , r G R and a G Z/dZ, satisfy the conservation rule 

d-\ d-l 
n=0 m=0 

where both sums do not depend on r and a. 
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The matrix F ra can be factorized as 

F ra = D ra F, F := F 00 (92) 
where D ra is the d x d diagonal matrix with the matrix elements 

(Y) \ — m(d-m)a/2+((i-i) 2 r/4-m((l-l)r/2c /gox 
V- L -'ra/mn • 1 u m,n \ J - } ) 

For fixed d, there are one (/-multiple infinity of Gaussian matrices D ra (and thus F ra ) distinguished by 
a E Z/gB and ret. On the other hand, F is the well-known ordinary DFT matrix. The matrix F was 
the object of a great number of studies. The main properties of the ordinary DFT matrix F are summed 
up in the appendix. 

6.2.2. Hadamard Matrices 



The matrix F ra defined by (85) is unitary. The modulus of each of its matrix elements is equal to 
1 / \/~d. Thus, F ra can be considered as a generalized Hadamard matrix (we adopt here the normalization 
of Hadamard matrices generally used in quantum information and quantum computing) [26-31]. 

In the case where d is a prime number, Equation (88) shows that the matrix (F ra )^F r b is another 
Hadamard matrix. However, it should be mentioned that, given two Hadamard matrices M and N, the 
product M^N is not in general a Hadamard matrix. 

6.2.3. Trace Relations 



The trace of F ra reads 

trF ra = e^ d ^ 2 ^-^S(u,v,w) 
yd 

where S(u, v, w) is given by (78) with 

u := 2 — a, v := d(a — r) + r, w := d 



(94) 



(95) 



Note that the case a = 2 deserves a special attention. In this case, the quadratic character of tr F r 
disappears. In addition, if r = we get 



tr F 02 = Vd 



(96) 



as can be seen from direct calculation. 



6.2.4. Diagonalization 

It is a simple matter of calculation to prove that 



(F^) 1 " V m F ra = q 



(d-l)(r+a)/2 



(q 1 
q 2 







(97) 



V ... q d J 
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where the matrix 








q a 
q 2a 





V £>(<*-!)'- o o 



q 



o 



(d-l)a 

o y 



(98) 



represents the linear operator v ra defined by (68). Therefore, the matrix F ra reduces the endomorphism 
associated with the operator v ra . 

Concerning (97) and (98), it is important to note the following conventions. According to the tradition 
in quantum mechanics and quantum information, the matrix V ra of the operator v ra is set up on the basis 
B 2 j + \ ordered from left to right and from top to bottom in the range = \d — 1), — 1) = \d — 
2), . . . , \j, —j) = |0). For the sake of compatibility, we adopt a similar convention for the other matrices 
under consideration. Thus, the lines and columns of F ra are arranged in the order d — 1, d — 2, . . . , 0. 

6.2.5. Link with the Cyclic Group 

There exists an interesting connection between the matrix X and the cyclic group Cd [24-26]. Let us 
call R a rotation of 2-k /d around an arbitrary axis, the generator of Cd- Then, the application 



C d ^ {X n :n = 0,l,...,d-l} : R ^ X 



(99) 



defines a rf-dimensional matrix representation of Cd- This representation is the regular representation of 
C d - Thus, the reduction of the representation {X n : n — 0, 1, . . . , d — 1} contains once and only once 
each (one-dimensional) irreducible representation of C d . 

6.2.6. Decomposition 



The matrix V ra can be decomposed as 



P r XZ a 



(100) 



where 



/l 



1 






1 



and 



/0 1 




X :-- 




1 




\1 



\0 













(101) 



1 

0/ 



Z := 



r(d-l)r J 

/l 
q 

q 2 

Vo 



\ 






-,d-l 



(102) 
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The matrices P r , X and Z (and thus V ra ) are unitary. They satisfy 

V ra Z = gZV ra (103) 
V ra X = f a XV ra (104) 

Equation (103) can be iterated to give the useful relation 

(V ra ) m Z n = q mn Z n (V ra ) m (105) 
where m,n G Z/dZ. Furthermore, we have the quasi-nilpotency relations 

e- i7r(d - 1)r (V r0 ) d = Z d = I d (106) 
(the relations (106) are true nilpotency relations when r = 0). More generally, we obtain 

Vn G Z/dZ : (V ra ) n = g -™(«-iW2( Vr0 ) n Z an (V ra ) d = <>(*-!)('+«») i d (1Q7) 
in agreement with the obtained eigenvalues for V ra (see Equation (71)). 
6.2.7. Weyl Pairs 

For r = a = 0, Equations (103) and (106) show that the unitary matrices X and Z satisfy the 
g-commutation relation 

XZ = gZX (108) 

and the nilpotency relations 

X d = Z d = I d (109) 

Therefore, X and Z constitute a Weyl pair (X, Z). Note that the Weyl pair (X, Z) can be defined from 
the matrix V ra only since 

X = V 00 , Z = (V o) f Voi (110) 

which emphasize the important role played by the matrix V ra . Note also that according to (97), we have 

F f XF = gZ (111) 

that proves that X and Z are related by the DFT matrix. 

Weyl pairs were introduced at the beginning of quantum mechanics [38] and used for building 
operator unitary bases [39] . The pair (X, Z) plays an important role in quantum information and quantum 
computing. In these fields, the linear operators corresponding to X and Z are known as flip or shift and 
clock operators, respectively. For d arbitrary, they are at the root of the Pauli group, a finite subgroup of 
order d 3 of the group U(d) for d even and SU(d) for d odd [30,31]. The Pauli group is of considerable 
importance for describing quantum errors and quantum fault tolerance in quantum computation (see 
[40-43] and references therein for recent geometrical approaches to the Pauli group). The Weyl pair 
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(X, Z) turns out to be an integrity basis for generating the set {X a Z 6 : a, b G Z/dZ}. The latter set 
constitutes a basis for the Lie algebra of the unitary group U(d) with respect to the commutator law. This 
set consists of d 2 generalized Pauli matrices in d dimensions [30,31]- In this respect, note that for d = 2 
we have 

X = cr x , Z = a z , XZ = -ur y , X°Z° = a (112) 
in terms of the ordinary Pauli matrices a = I 2 , cr x , cr y , and cr z . 
6.2.8. Link with a Lie Algebra 

Equation (105) can be particularized to give 

X m Z n = q mn Z n X m , (m,n)eN 2 (113) 

Let us define the operator 

T( ni ,n 2 ) := q^Z n ^X n \ (n 1; n 2 ) G N 2 (114) 
It is convenient to use the abbreviation 

(n 1 ,n 2 ) = n T n = T {nuri2) (115) 
The product T n T m is easily obtained to be 

T m T n = q-? mxn T m+n (116) 

where 

m x n := min? — m 2 ni, m + n = (mi + ni,m 2 + n 2 ) (117) 
The commutator [T m , T n ], 

[T m , T n ] = -2i sin (j-m x n) T m+n (118) 

follows at once from (116). The operators T m can be thus formally viewed as the generators of the 
infinite-dimensional Lie algebra (or sine algebra) investigated in [44,45]. 

7. Closing Remarks 

We used the representation theory of the symmetry groups SU(2) and SU(1,1) to describe the 
generalized oscillator algebra A K and the two phase operators E d and introduced in [16]. The phase 
eigenstates of E d and E^ were thus understood in terms of finite-dimensional and infinite-dimensional 
representations of SU(2) and SU(1,1), respectively. In the case of those representations of SU(2) for 
which the dimension is a prime integer, our approach led us to derive MUBs as eigenbases of the phase 
operator E d (with d prime), opening a way for further results on unitary phase operators associated with 
Lie groups. 
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The unitary phase operator E d defined via 



(119) 



leads to a polar decomposition of the algebra su(2) in the scheme {J 2 , E d }, which is an alternative to 
the familiar quantization scheme { J 2 , J3} of angular momentum theory. The { J 2 , E d } scheme and the 
{J 2 , v ra } scheme of [22-32] are related by (74). In the case of the noncompact Lie algebra su(l,l), the 
phase operator E^ is non-unitary and given by 



Although this does not correspond to a true polar decomposition (because E^ is not unitary), it yields a 
scheme {K 2 , E^}, which is an alternative to the canonical scheme {K 2 , K 3 } developed for su(l,l) by 
Bargmann and most of other authors. We hope to further study this new scheme from the point of view 
of the representation theory and the Wigner-Racah algebra of SU(1,1). 

As far as the applications of the new SU(2) and SU(1,1) phase states derived in Section 3 are 
concerned, let us mention that, besides the two applications (to mutually unbiased bases in section 5 
and to discrete Fourier transform in Section 6) discussed in our paper, we can mention other potential 
applications. Our phase states can be useful for various dynamical systems (e.g., the Morse system for 
the SU(2) states as well as the Poschl-Teller system and the repulsive oscillator system for the SU(1,1) 
states). We can also mention a possible application of the quadratic discrete Fourier transform to discrete 
linear canonical transforms and to Hadamard matrices in connection with the production of geometric 
optics setups. Some of these further potential applications are presently under consideration. 
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Appendix: Properties of the Ordinary DFT Matrix 

The ordinary DFT — also called the finite Fourier transform — is the linear transformation of the 
complex d-dimensional Hilbert space C d onto itself, that is represented by the matrix F whose elements 
are given by 

(F) mn := -±=q mn = ^f^ mnld (121) 

with m,n = 0, 1, . . . , d — 1. The elements (F) mn are periodic in m and n modulo d (so that F can 
be stitched into a torus), but we shall consider the fundamental interval to be < m,n < d — 1. 
Equation (121) follows from (85) and (92). The matrix F corresponds to the transformation (89) with 

1 d-l 1 d-l 

* m=0 * n=0 

Note that in the physics literature it is more common to find the definition (121) with a minus sign in the 
exponent; of course, the results obtained with the two conventions are equivalent. 

The Fourier matrix F has several well-known properties. It is symmetric and unitary. In addition 
it satisfies 



(123) 
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Because F is unitary, its eigenvalues must be on the unit circle S 1 , and since it is a fourth root of unity, 
so are its eigenvalues, to be denoted by 

<^:=i* = e- fe / 2 e{l,i,-l,-i} (124) 

for k — 0, 1, 2, 3. This divides the space C d into four Fourier-invariant, mutually orthogonal subspaces 
whose dimensions N iPk exhibit the modulo-4 multiplicities of the eigenvalues </?*,.. Of course, we have 

3 3 3 

d = Y< N ^ t rF = J> fcA ^' detF = J](^ fc )^ (125) 

k=0 k=0 k=0 

For d = 4 J + k with k = 0, 1,2,3 and J E N, the multiplicities, traces and determinants of the 
submatrices of F associated with each eigenvalue are given by: 

dimension multiplicities 

d= V 9 o=l V 9 i = i f2= — 1 y?3= — i tr F det F 

4J J+l J J J-l 1 + i i(-l) J ( 

4J+1 J+l J J J 1 (-1) J 

4J+2 J+l J J+l J -(-1) J 

4J+3 J+l J+l J+l J i -i(-l) J 

(see for example [46] noting that the DFT matrix there is the complex conjugate of the DFT matrix here). 

Since N Vk ps d/A, there is wide freedom in choosing eigenvector bases within each eigenspace. 
Finding a "good" eigenbasis is of interest to define fractional powers of the DFT matrices, which 
constitute the abelian group of elements {F^}, for real v modulo 4 [47], that would contract, for d — > oo, 
to the fractional Fourier integral transform. The fractionalization of the Fourier integral transform was 
defined in 1937 by Condon [48] at the suggestion of von Neumann, rediscovered in other contexts 
[49,50], and is currently of importance for signal analysis and image processing through the fast Fourier 
transform algorithm [51-54]. The integral kernel of the fractional Fourier integral transform can be 
expressed as a bilinear generating function for Hermite-Gauss functions [55], 

vp n (x) := 1 e- x2 ' 2 H n {x) {111) 
y / 2 n n\ tJtx 

where H n (x) are the Hermite polynomials of degree n E N in x, which are the eigenfunctions of the 
Fourier integral transform J 7 , 

^ i r+oo 

tt n (x) = 0F*„)(aO = -=r- / dx'e ixx ' * n (x') = i n V n (x) (128) 

V 271 " J-oo 

The integral kernel of the fractional Fourier integral transform [47] is then obtained as 

„, 1 / xx' — Ux 2 +x' 2 ) cos \txu\ 
F v (x,x') := exp i 2 \ , > 2 —) (129) 

^Trsinivrz/ V sm 2^ 7 

oo 

= J2^)^ /2 ^n(x') (130) 
n=0 
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for < v < 2, with the limits 

F°(x,x') = 5(x-x') 1 F 2 (x,x') = 5(x + x') (131) 
These kernels are unitary, 



F- v (x,a/) = F v (x',x) (132) 

and form a one-parameter group 

/+oo 
dx' F Ul (x,x')F U2 (x',x") = F Ul+U2 {x,x") (133) 
•00 

with v modulo 4. 

To fractionalize the DFT matrix F, one will be naturally interested finding <i-point functions $n (m) 
that are "good" discrete counterparts for the Hermite-Gauss functions ^> n (x) in (127); in particular that 
they be analytic and periodic functions of m. Mehta [56] has proposed the following functions: 

$£°(m):= ex P (-^ + m) 2 )# n (J-^ + m)) (134) 

e=-oo 

that we call Mehta functions. These have the desired properties and 

F$W=i n $W (135) 

for n6N and where is the column vector of components (m) . Of course, there cannot be more 
than d linearly independent vectors in C d , so we may take the subset {3>i^ : n — 0, 1, . . . , d — 1}. Prima 
facie, it is not clear whether this subset is linearly independent and orthogonal, or not - Mehta [56] left 
unresolved their orthogonality, which was lately described thoroughly by Ruzzi [57]. The departure from 
strict orthogonality of the vectors of the Mehta basis was investigated in [58]; the departure is small for 
low values of n and gradually worsens up to d — 1 . 

Indeed, there is wide freedom in choosing bases for C d when the sole requirement is that they be 
eigenbases of F, satisfying (135). Labelling these eigenvectors by their four Fourier eigenvalues Lfk, 
and within each of these eigenspaces by j — 0, 1, . . . , N Vk — 1, we denote them by {T^ ,j)(m)}, 
periodic in m modulo d [46,58]; and we assume that they are complete in C d and thus have a dual basis 
{^(fk,j)( m )} periodic in m, such that 

d-l 

Yl %k,d m ) T (w,i')( m ) = s k,k>Sjj> (136) 

m=0 

and 

N Vk -l 

Yl T fej)H^j)( m ') = ( n ft)mm' (137) 
3=0 

where 11^ is the projector matrix on the Fourier subspace tp k . Associated with this basis {T}, one may 
define the corresponding 'T-fractionalized DFT matrices' F^. with elements 

3 N ¥k -1 

(ESrW := E E T (^,i)(^) ^ mk)v/2 % k d™') d38) 

fc=0 j=0 
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where we use the compound index n = 4j ' + k to order the vectors, as if it were the 'energy' label in 
the Mehta functions (134). In this way, the vectors of the {T} basis are eigenvectors of a number matrix 
N-f with elements 

3 N Vk -l 

(NtW : = E E T (^)( m ) ^ +k ) % k ,d m ') (139) 

fc=0 j=0 

In other words 

NrT w =nT w (140) 
The matrix N-f has the virtue of being the generator of the T-fractional Fourier matrices, 

F£ = exp(iivrz/N T ) (141) 

for v modulo 4. 
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